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Theoretical estimates for the half life of neutrinoless double beta decay (Ovfip) in candidate nuclei 
are affected by both particle and nuclear physics uncertainties, which may complicate the interpre¬ 
tation of decay signals or limits. We study such uncertainties and their degeneracies in the following 
context: three Ov/3/3 nuclei of great interest for large-scale experiments ( 76 Ge, 130 Te, 136 Xe), two rep¬ 
resentative particle physics mechanisms (light and heavy Majorana neutrino exchange), and a large 
set of nuclear matrix elements (NME), computed within the quasiparticle random phase approxima¬ 
tion (QRPA). It turns out that the main theoretical uncertainties, associated with the effective axial 
coupling gA and with the nucleon-nucleon potential, can be parametrized in terms of NME rescaling 
factors, up to small residuals. From this parametrization, the following QRPA features emerge: 

(1) the NME dependence on gA is milder than quadratic; (2) in each of the two mechanisms, the 
relevant lepton number violating parameter is largely degenerate with the NME rescaling factors; 
and (3) the light and heavy neutrino exchange mechanisms are basically degenerate in the above 
three nuclei. We comment on the challenging theoretical and experimental improvements required 
to reduce such particle and nuclear physics uncertainties and their degeneracies. 


I. INTRODUCTION 

The search for the neutrinoless mode of double beta decay (Ozz/3/3) in different (Z , A) candidate nuclei, 

( Z , A ) —^ (Z + 2, A ) + 2e , (1) 

represents a major research program in experimental neutrino physics M- From the theoretical viewpoint, this 
decay mode provides a unique probe of the Dirac or Majorana nature of neutrinos [6|] and, in general, of lepton 
number violation (LNV) processes 0-H3- Indeed, the decay may be mediated not only by Majorana neutrinos with 
light (sub-eV) masses, but also by other LNV mechanisms involving new particle physics at higher mass or energy 
scales 

For a single LNV mechanism (labeled by the index j) occurring in a candidate nucleus (labeled by the index i), the 
Oi//3/3 decay half life T) can be expressed as 

Tr'=Gi\Mi\\M? , (2) 

where G\ is a kinematical phase-space factor, M\ is the nuclear matrix element (NME) encoding the nuclear dynamics 
of the process, and X 3 is a LNV parameter encoding the particle physics aspects of decay. In principle, by means of 
independent Ozz/3/3 decay observations in different nuclei, one may hope to disentangle the underlying particle physics 
mechanism and the associated A 3 value [TTHIi]. Unfortunately, the relatively large nuclear model uncertainties 
affecting M 3 [Tfl [Tot | make this program quite difficult in practice. 

Several studies have investigated the stringent conditions on the M 3 uncertainties, under which various particle 
physics mechanisms may—or may not—be disentangled with future, multi-isotope Oiz/3/3 decay data (see [7i. !TT1 ml Sol- 
1251 for an incomplete bibliography). A new twist in this field has been provided by recent discussions on the axial 
coupling gA, which could be significantly suppressed with respect to the vacuum value < 7 ™ c cs 1.27, due to nuclear 
medium and other “quenching” effects. In particular, qa values in the reference range often used in the past literature 


1 <g A < 1.27 , (3) 

are being increasingly questioned by phenomenological studies. In fact, in the framework of the quasiparticle random 
phase approximation (QRPA), best-fit values g A < 1 were obtained in [28] by a joint analysis of the experimental 
data available for the electron capture (EC), single beta (/3) and two-neutrino double beta (2i/2/3) processes in the 
100 Mo and 116 Cd nuclei. 
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For the same weak processes and nuclei, such early indications for gA < 1 l28j were later confirmed in another 
QRPA approach [53. Recently, a marked preference for values of gA below unity has also been found in a different 
approach based on the interacting boson model (IBM) j30l43^ |. where a possible dependence of gA on the atomic 
number has been noted [3(|. See also [3^ for a recent discussion of quenching variations in the chiral effective field 
theory approach. In the absence of a deep theoretical understanding of quenching effects UM , these findings suggest 
that the usual reference range in Eq. ([3| should be conservatively extended somewhat below unity, e.g., in the range 
0.8 < g A < 1.27 0. 

In general, small (or uncertain) values of gA may strongly affect the half-life estimates via Tf 1 oc |Mf| 2 oc g\, thus 
making even more difficult to constrain the particle physics mechanism and its LNV parameter [351. 136j|. However, 
in the adopted QRPA framework (0, the (^-dependence of \M?\ is known to be milder than quadratic, since gA 
variations may be partly traded by shifts of another free parameter—the particle-particle strength g pp [38[—via a 
joint fit to reference 2vfiP data; see, e.g., [53, HD, 0, 38, 0. Therefore, it makes sense to revisit the problem of 
determining the particle physics mechanism and its LNV parameter, by allowing gA < 1 within the QRPA. 

To this purpose, we consider in the following three nuclei ( 76 Ge, 130 Te, 136 Xe) and two 0v/3[3 mechanisms mediated 
by light (L) and heavy ( H ) Majorana neutrino exchange, within an updated QRPA approach. We show that, even if 
the decay half lives were accurately measured, the current nuclear model uncertainties (mainly related to gA and to 
the nucleon-nucleon potential) would lead to a degeneracy between the LNV parameter and the NME errors in each 
mechanism, as well as between the mechanisms themselves. Although limited to a few representative nuclei and Oi//3/3 
decay processes, these results highlight the severe conditions and the challenging improvements needed to (partially) 
lift such degeneracies in the future. 

Our work is structured as follows. In Sec. II we introduce the notation and conventions for the two (L and H) decay 
mechanisms. In Sec. Ill we discuss the QRPA calculation of the nuclear matrix elements and the parametrization of 
the associated uncertainties. In Sec. IV we perform a statistical analysis of prospective data, showing the degeneracy 
of particle and nuclear physics uncertainties, both within each mechanism and between the two mechanisms. We 
briefly summarize our results in Sec. V. 


II. LIGHT AND HEAVY NEUTRINO EXCHANGE: NOTATION AND CONVENTIONS 


In the following, we shall study Ozz/3/3 decay in three representative nuclei of great interest for large-mass projects, 

i = 1, 2, 3 = 76 Ge, 130 Te, 136 Xe , (4) 

and two reference LNV mechanisms, mediated by either light ( L ) or heavy (H) Majorana neutrino exchange [40l [4lj| . 


j = 1, 2 = L, H . (5) 

We refer the reader to [M3 for recent discussions of the particle physics dynamics of the L and H mechanisms, and 
to (42] for an approach interpolating between these two cases. 

Here we just remind that, for the L mechanism, the LNV parameter can be expressed as [1]: 


X L = m,f 3/3 = 


E 

h =1 


\Ueh\ 2 e^ h m h 


( 6 ) 


where rrih are the masses of the three known light neutrinos U e h are their mixing matrix elements with v e , and <f>h 
are unconstrained Majorana phases (one of which can be rotated away). For the H mechanism, the LNV parameter 
can be expressed as (see, e.g., 0): 


X H = Mgp = m e 


E Vekfj** 

k>4 


1Tl p 

M k 


(7) 


where M k are the masses of possible heavy neutrinos beyond the known ones (assuming M k m p ), U ek are their 
mixing matrix elements with is e , and are further Majorana phases. 

The L and H mechanisms are characterized by the same phase space (4lj | . 

Gf = Gf = G t , (8) 

which, in the conventions of [37], embeds a factor l/mjji, so that its units are [Gi\ = y -1 eV -2 , while [A- 7 ] = eV. 
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TABLE I: Phase space values [43j] and 90% C.L. half-life limits for the three nuclei considered in this work. See the text for 
details. 


i 

Gi ( y - 1 eV 2 ) 

7* 

Ti (y) 

Experiments 

76 Ge 

2.21 x 10 -26 

25.656 

> 3.0 x 10 25 

GERDA + IGEX + HdM [44] 

130 Te 

1.33 x 10 -25 

24.876 

> 4.0 x 10 24 

CUORE-O + CUORICINO [47] 

136 Xe 

1.36 x 1CT 25 

24.865 

> 3.4 x 10 25 

KamLAND-Zen + EXO-200 [49] 

Finally, 

we linearize Eq. © by taking logarithms as in [b| 53,53 > 




n = 7* - 

- H - V , 

(9) 

where 


O 

bio 

O 

II 

(Ti/ y) , 

(10) 



I 

II 

0 

CTQ 

0 

[Gi/(y _1 eV -2 )] , 

(11) 



vl = log 10 

wi\ , 

(12) 



^ = log 10 

(A-’/eV) . 

(13) 


Table I reports the numerical values of Gi (and of 7 j) as taken from (43[, together with the most stringent experi¬ 
mental lower limits for the half lives Xi as quoted by GERDA (44| for 76 Ge (in combination with IGEX [45| and HdM 
0), by CUORE 47 1 for 130 Te (in combination with CUORICINO [HJ]), and by KamLAND-Zen [4£j for 136 Xe (in 
combination with EXO-200 0 ). Note that the G, values of (43] have been rescaled by the fourth power of 1.27, for 
consistency with our conventions [see Eq. (TTH) below]. 

We emphasize that, although the L and H mechanisms share the same phase-space factor G;, their dynamics is 
quite different. The exchange potentials behave like 1/r and delta functions in the L and H cases, respectively, 
inducing significant differences in the multipole decomposition of the corresponding matrix elements, as well as in 
in their sensitivity to gA changes (not shown). Furthermore, the H mechanism is more sensitive than the L one to 
the choice of the nucleon-nucleon potential (see Sec. Ill B). Given these intrinsic dynamical differences, one does not 
expect a priori that the L and H mechanisms are phenomenologically degenerate, as they turn out to be. 


III. QRPA NUCLEAR MATRIX ELEMENTS AND UNCERTAINTIES 
A. Nuclear Matrix Elements 

The 0v/3[3 nuclear matrix element M = M J - consists of the Fermi (F), Gamow-Teller (GT) and tensor (T) parts 
which, in the adopted conventions, read [371 kill.1521 : 

M =(rl) 2 (-(^ +MCT -"f) ■ (14) 

Note that the effective axial coupling gA enters not only in the prefactor and in the Fermi matrix element, but also in 
the calculation of the GT and tensor constituents, due to a consideration of the nucleon weak-magnetism terms fill ]. 

In this work, six sets of nuclear matrix elements M? have been calculated for each of the three nuclei in Eq. GD and 
of the two mechanisms in Eq. ([5]). In particular, for both light and heavy neutrino mass mechanisms, the calculation is 
based on the QRPA with partial restoration of isospin symmetry (53$. The nuclear radius R = roA 1 / 3 with ro = 1.2 fm 
is used. For each pair (i, j), the NME set includes 18 variants, according to three different sizes of the single-particle 
space (small, intermediate, large), two different types of nucleon-nucleon interaction (charge-dependent (CD) Bonn 
and Argonne) [52], and three different values of gA , 

g A = 1.27, 1.00, 0.80 , (15) 

which are representative of unquenched, quenched, and strongly quenched axial couplings. We remind that, for 
each calculation, the g pp value is fixed by imposing that the corresponding (theoretical) 2v(3(3 half-life equals the 
experimental one in each nucleus [37l.l5ll]. Summarizing, the following analysis is based on a total of 108 Oiz/3/3 matrix 
elements, calculated within a QRPA framework which reproduces three 2v(3(3 half-lifes by construction. 
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FIG. 1: Scatter plot (in logarithmic scale) of the nuclear matrix elements \M?\, for each pair of the three 0 vj3fi candidate 
nuclei, and for both heavy (H) and light (L) Majorana neutrino exchange. In the latter case, the NME have been multiplied 
by a factor 100. See the text for details. 


Figure 1 shows the scatter plots of these 108 NME values in logarithmic scale, which visualize the correlations 
between pairs of the rjj parameters in Eq. m- NME variants are distinguished by different marker types: blue and 
red for L and H mechanisms; full and hollow for Argonne and CD-Bonn potentials; and squares, circles and triangles 
for gA = 0.80, 1.00 and 1.27, respectively. Single-particle space sizes are not explicitly distinguished (see also below). 
Note that, for graphical convenience, the Mf values for the L mechanism have been multiplied by a factor 100. This 
figure reveals a very strong correlation of the NME’s, which is especially evident in the plane charted by 136 Xe and 
130 Te. In practice, it appears that the theoretical uncertainties associated to NME variants can be largely absorbed 
by overall scaling factors. 

The strong correlations emerging in Fig. 1 imply two different types of degeneracies between nuclear and particle 
physics aspects of the decay, at least in the (i, j) sets considered herein: (a) for a given decay mechanisms, the nuclear 
model uncertainties are degenerate with the LNV parameter A J 17]; and (b) the two different mechanisms (L and 
H) are largely degenerate with one another [2(j. With respect to [13, H3], we sharpen these statements by using a 
convenient parameterization and statistical treatment of the NME uncertainties within the QRPA (including cases 
with gA < 1), as discussed below. 

A final remark is in order. Variations of the single-particle space size (from small to intermediate and large size) 
produce relatively small NME changes, partly orthogonal to the main degeneracy directions of Fig. 1. These changes 
do not reveal a specific pattern: e.g., it turns out that, for 76 Ge, the NMEs are slightly low for intermediate space 
size (with respect to small or large sizes), while the opposite happens for 130 Te and 136 Xe; for the latter two nuclei, 
the NMEs values for small and large sizes somewhat differ, which they are rather close for ,6 Ge. In the absence of a 
compelling pattern emerging from different single-particle space sizes, we omit their distinction in Fig. 1, and adopt 
the conservative viewpoint that the associated NME variations are basically uncorrelated among the three nuclei. 
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TABLE II: Numerical values for the empirical parametrization of nuclear model uncertainties in Eq. EH, for the two L and H 
mechanisms. 


3 

g( 

V J 2 

V 3 

a? 

0 j 

cr- 7 

L 

0.600 

0.504 

0.267 

0.458 

0.021 

0.032 

H 

2.400 

2.364 

2.135 

0.544 

0.089 

0.025 


B. Parametrization of uncertainties 

We find that the spread of the numerical 77 3 values can be fully covered by the following empirical parameterization, 
which incorporates both the strong linear correlation among the rj 3 and the residual transverse scatter visible in Fig. 1: 

vj = fi + a 3 ( g A - 1 ) + s0 3 ± a 3 , (16) 

where s = +1 (—1) for the CD-Bonn (Argonne) potential, and the parameters fjj, a 3 , 0 3 and a 3 are given in Table II. 
In the above equation, fj 3 represents a sort of “central value” for the set of r] 3 values, while a 3 (gA — 1) + s0 3 represents 
the systematic theoretical bias due to admissible variations of gA with respect to the unit value, and of the nucleon- 
nucleon potential. In terms of NMEs (\M\ = lCf 7 ), the bias acts as an overall (z-independent) NME rescaling factor 
for the three nuclei. Finally, ±a 3 represents the residual range which is not covered by the previous bias, including 
variations due to the basis size (small, intermediate, large). [Actually, the a 3 values covering the r/ 3 spread depend 
slightly on the index i\ we neglect these tiny variations, and conservatively take the largest value for a 3 .} 

The above parametrization indicates that, within the QRPA, the functional dependence of the Ov00 NME on gA 
is significantly milder than the naive quadratic expectations (| M 3 | oc g ^), as already noticed in ( 2 ^, |3^ [3^]. We recall 
that, within the QRPA approach, the g pp parameter is adjusted to fit the 2v00 decay rate, and that both the Ou00 
and 2v00 NME (M 0v and M 2v ) decrease with decreasing gA or with increasing g pp . Then, if gA decreases, the g pp 
parameter must decrease as well, in order to keep M 2u at the value fixed by the 2i/00 half life T 2l/ , as shown in Fig. 2. 
As a consequence, also the change in the matrix element M 0v is smaller than one might at first expect (see, e.g., 
[H|). In particular, Eq. (fl6l) suggests that, for Ou00 decay, the effective NME dependence on gA is close to be linear 
(| M 3 1 oc gA) rather than quadratic, at least for relatively small values of the difference gA — 1, and within the rough 
approximation a 3 ~ 1/2 ~ ln(10). In this sense, the impact of the large gA uncertainties in the interpretation of Qv00 
data [13, HI, Hlj may be effectively reduced (although not eliminated) within the QRPA approach [H, HI, HI'] • 



FIG. 2: The relation between the weak axial-vector coupling parameter gA and the g pp parameter, determined from the 
measured 2v00 half-life of 76 Ge, 130 Te and 136 Xe. The results refer to the case with Argonne potential and small size for the 
single-particle model space. Similar results hold for other choices of the potential or model space (not shown). 
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We remark that, for the sake of simplicity, gA has been assumed to be the same in the three considered nuclei. This 
seems to be an acceptable starting point for a phenomenological analysis, since the origin and amount of quenching 
are not well known. For instance, quenching effects might be assigned to the A-isobar admixture in the nuclear 
wave function, or to the shift of the GT strength to higher excitation energies due to short-range tensor correlations. 
Model-space truncation can also exclude strength that may be pushed to high energies, and the omission of two- 
body currents can also leave excitations unaccounted for. Moreover, quenching effects might be different for different 
multipoles and, if associated with exchange currents, might be smaller for light nuclei. In the absence of a clear picture 
for these effects, we have assumed that the same value of gA (either 0.8, 1.0, or 1.27) applies to all multipoles in the 
three medium-heavy nuclei consider herein. Of course, this simplified assumption may be revisited in future and more 
refined analyses of the quenching phenomenon within the Or^/3/3 decay context. In perspective, one should build a 
general theory of quenching in the nuclear medium, and constrain systematically the theory with data from different 
nuclear processes linked to 0u/3f3 decay, so as to reduce the effects of the gA uncertainty (see, e.g., the discussion in 
[l^lU). At present, however, we must accept theoretical uncertainties (at least) as large as in Eq. (flbl) . 

So far, we have mainly discussed the sensitivity to gA variations, characterized by the a? parameters. Let us now 
comment on the other parameters ft and ft . The value of ft , which characterizes the NME sensitivity to the choice 
of the nucleon-nucleon potential, turns out to be much larger in the H mechanism than in the L one (by a factor 
of about four), as anticipated at the end of Sec. II. Following the remarks at the end of Sec. Ill A, the theoretical 
uncertainty ±ft is treated as a “one-standard-deviation range”, covering those residual uncertainties (including the 
single-particle space ones) which are not included in the “bias” term a? (gA — 1) + sft . This definition is conservative, 
because it allows to cover [via Eq. (fTBl) ] all the NME in Fig. 1, and not only 68% of them. 


C. Degeneracies in terms of observable half lives 

The results of Sec. Ill B allow to visualize the degeneracies mentioned at the end of Sec. Ill A in terms of observable 
quantities, i.e., the Qvfifi decay half lives Xj expected in different nuclei and for different underlying mechanisms. Since 
the Ti are usually represented in logarithmic scale, we shall base our discussion directly on Ti = log 10 (Tj/y). From 
Eqs. © and (usd, the Ti can be expressed as 

Ti = 7 * - 2 77 ^ - 2ft ± 2ft , (17) 

where we have defined a “rescaled” LNV parameter ft as, 

ft = ft + ft (gA — 1 ) + sft . (18) 

The above equations clearly show the degeneracy between the particle physics parameter ft and the systematic QRPA 
uncertainties parameterized by a?(gA — 1) + sft: variations of the latter term can be traded for opposite changes in 
the LNV parameter, without affecting the observable Tj. 

Figure 3 shows the theoretical expectations in the planes charted by pairs of half lives (t*,, t/j) in different nuclei, 
together with the current experimental lower limits as reported in Table I. The half lives are correlated via Eq. m , 
which implies 

Tk - T h = - 7 ft - 2 (fg - rf h ) (19) 

for j = H, L , up to residual errors (±2cr J ), shown as crosses in Fig. 1. The position of each cross is irrelevant: the 
associated errors are the same at any point on the slanted lines—which should thus be thought as “error bands”. This 
figure illustrates the second kind of degeneracy mentioned at the end of Sec. Ill B, namely, the near indistinguishability 
of the L and H mechanisms, which is especially evident in the rightmost panel, where the half-life expectations for 
the L and H cases are almost coincident. In the other two panels on the left, the L and H lines are visually separated, 
but they can overlap within error bars. As far as the separation between the slanted lines in Fig. 3 remains smaller 
or comparable than the theoretical error bars, future experimental data on the half lives (no matter how accurate) 
will not be able to tell the L from the H mechanism. In the following section, we shall quantify the two kinds of 
degeneracies by performing an analysis of prospective data. 


IV. STATISTICAL ANALYSIS OF PROSPECTIVE Ou/3/3 DATA 


There are good prospects to improve the current limits (Table I) on the Oi ft(3 half-life in the candidate nuclei ( 76 Ge, 
130 Te, 136 Xe). In particular, in the next decade, upgrades of the current experiments or new planned projects should 
be able to explore an additional order of magnitude in Ti [ZHlI and hopefully find evidence for the decay. Below, we 
shall optimistically assume that a 0^/3/3 decay signal is found in each of the three nuclei. 
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FIG. 3: Correlation plot of expected half lives in different pairs of nuclei (slanted lines) together with residual theoretical errors 
(crosses) for the L and H decay mechanisms (blue and red lines, respectively). Also shown are the current experimental limits, 
as reported in Table I. See the text for details. 


A. Reference scenario 

We assume that the measured ' 6 Ge half life takes a reference value of 10 26 y. The half lives for 130 Te and 136 Xe in 
the L and H mechanisms are then obtained via Eq. (Hhd. We also assume that each X) is measured within ±20% at 
ler, corresponding to ±0.08 accuracy on n. Our prospective data sets are thus given by the following central values 
and errors {ff ± sj) for the L and H cases, respectively: 

{ if ±s[ = 26.000 ± 0.080 , 

if ± sf = 25.412 ± 0.080 , (20) 

if ± sf = 25.875 ± 0.080 , 

i ff ± sf = 26.000 ± 0.080 , 

ff ± sf = 25.292 ± 0.080 , (21) 

ff ± sf = 25.739 ±0.080 . 

Figure 4 shows the above data sets (black dots with crossed error bars) overlaid on the same theoretical predictions 
(colored slanted lines with errors) for the L and H cases as in Fig. 3. We discuss below the implications of these data 
on the degeneracy issues, by means of a statistical analysis. 








FIG. 4: As in Fig. 3, but with lower limits replaced by prospective determinations of the half lives (black dots and crosses). 
See the text for details. 


B. Degeneracy between LNV parameter and QRPA uncertainties for fixed Ov/3/3 decay mechanism 


Let us assume one of the two Majorana neutrino exchange mechanisms (either L or H) as the “true” one for the 
Ovpfj decay in the three nuclei. Given the previous QRPA calculations, the prospective data sets, and their associated 
uncertainties, we aim at determining the associated LNV parameter (either /i L or g H ) and its errors. 

To this purpose, we consider the following y 2 function in terms of the “rescaled” parameter £l 3 in Eq. 




xV) = ^ 


in-n) 


J\2 


f=l ( s i) 2 + (2ad) 2 


3 

E 

i =1 


( 7 i - 27 ^ - 2 fX 3 - Tj) 2 

{si) 2 + M 2 


( 22 ) 


Minimization of y 2 provides the central value jjp c (at y 2 


0, by construction) and its error 5 J (at Ay 2 


1). We get: 


-0.772 ± 0.029 (L) , 
-2.572 ± 0.027 ( H) . 


The LNV parameter /V is then obtained from Eq. (18) as 

/i J = £l 3 c — oP{gA — 1) — sft ± (P . 


(23) 


(24) 


The above expression for g, 3 provides a useful breakdown of its uncertainties: from left to right, the terms following the 
central value J1 3 C represent, respectively, the systematic error due to gA variations from the unit value, the systematic 
bias due to the choice of the nucleon-nucleon potential (CD-Bonn vs Argonne), and the residual error from theory 
and data uncertainties. 
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FIG. 5: Upper panel: Light (L) Majorana neutrino exchange mechanism. The LNV parameters mpp (in eV), as derived from 
a fit to the prospective data in Eq. is shown as a function of qa (in the range qa G [0.8, 1-27]). The two curves refer to 
the CD-Bonn and Argonne choices for the nucleon-nucleon potential. The error bars attached to each line mark the size of the 
residual theoretical and experimental uncertainties from the fit. Lower panel: as above, but for heavy (H) Majorana neutrino 
exchange with LNV parameter Mpp, from a fit to the prospective data in Eq. m- 


Figure 5 represents the results of the above statistical analysis, in terms of the LNV parameters mpp and Mpp 
associated to the L and H mechanisms, as defined in Eqs. © and (|T|). respectively. The LNV parameters are shown 
(in logarithmic scale) as a function of gA , whose variation in the representative range gA £ [0.8, 1.27] induces the 
main source of uncertainty: the higher gA, the smaller mpp or Mpp. In decreasing order of relevance, the second 
source of uncertainty is represented by the nucleon-nucleon potential, whose twofold option (CD-Bonn vs Argonne) 
splits the LNV parameter estimates into two curves. Finally, the third and smallest source of uncertainties is induced 
by [see Eq. USD] and is shown as a representative error bar, on top of each curve. The total uncertainty affecting 
each LNV parameter is given by the “envelope” of all these errors which, for gA £ [0.80, 1.27], amounts to a factor of 
~ 2 for mpp, and to a factor of ~ 3 for Mpp (from minimum to maximum, in both cases). 

Overall uncertainties of a factor ~ 2 in the reconstructed value of mpp are definitely large, but not as dramatic as 
those of 0(10) discussed in [13,®]. This feature of the QRPA approach may thus be of interest for relatively robust 
estimates of the experimental sensitivity to mpp. We stress that, as mentioned in Sec. Ill B, in the QRPA the effective 
dependence of the NME on the gA parameter (M- ~ gA) is milder than in other frameworks (where M\ ~ g\), as a 
result of the stabilizing role of the u 2v/3/3 calibration” of the g pp parameter. 




10 


In any case, limiting the gA range in Fig. 5, by means of dedicated theoretical and experimental studies, will be 
a major step towards the reduction of the reconstructed LNV parameter uncertainties. One should also refine the 
understanding of the nucleon-nucleon potential, so as to bring the two splitted curves in Fig. 5 closer to each other. 
Such a long-term nuclear modeling program, although rather challenging, is warranted by the fundamental importance 
of the worldwide Oi '(3(3 decay search program. 


C. Degeneracy between light and heavy neutrino exchange mechanisms for Oi'/3/3 decay 

As already mentioned, the ensembles of expected half lives for the L and H mechanisms (slanted lines in Fig. 3) 
are very close to each other, as compared to current theoretical uncertainties and prospective experimental errors 
(see Fig. 4). This fact suggests that the two mechanisms are largely degenerate, i.e., they cannot be distinguished 
by data from the three candidate nuclei considered herein. One can quantify the degree of degeneracy as follows: 
the prospective data for one mechanism (say, L) are fitted with the predictions of the other mechanism (say, H), 
and viceversa. The value of xLin quantifies then the “degree of misfit”: the lower is Xmin, the more difficult is to 
distinguish the L and H mechanisms. By using a y 2 approach as in Eq. (1221) . we find that, as expected, the misfit 
is not statistically significant: Xmin — 1-1 i n both cases. Therefore, the two mechanisms are phenomenologically 
indistinguishable at the level of (Xmin )^ 2 ~ 1 a■ 

It should be noted that this small ler difference between the L and H scenarios is generated solely by the smallest 
sources of uncertainties (the data errors sj and the residual theoretical errors <r J ), while it does not depend on the value 
of gA or on the choice of the nucleon-nucleon potential (which can be absorbed by variations of the unknown LNV 
parameter Xj in the fit). Improving the latter two sources of uncertainties, despite being crucial to assess Xj , would 
not help to lift the L-H mechanism degeneracy. Breaking the degeneracy would require the challenging reduction of 
the non-parametric component of the theoretical error <t j (which is at the level of 10' Ti ~ 5-7% in our approach) and 
of the prospective experimental errors s* (which we have assumed to be at the ~ 20% level in this work). 

It is not obvious how the above stringent requirements can be achieved, even in the far future. In any case, we 
remind the reader that such conclusions refer to the specific decay mechanisms, candidate nuclei, and QRPA nuclear 
model considered in this work, and might thus be altered in a wider context. In general, favorable cases for the 
discrimination of any two decay mechanisms in a pair of candidate nuclei can be diagnosed, via correlation plots 
analogous to our Figs. 1 and 3, by the emergence of a significant “transverse” separation of the slanted error bands 
(see also gM). 


V. SUMMARY 

We have studied in detail a phenomenological scenario involving three candidate nuclei for 0vj3f3 decay ( 76 Ge, 
130 Te, 136 Xe), two representative particle physics mechanisms (light and heavy Majorana neutrino exchange, with 
LNV parameters mpp e Mpp), and a large set of nuclear matrix elements, computed within the quasiparticle random 
phase approximation. We have found that the main theoretical uncertainties, induced by the effective axial coupling 
gA and with the nucleon-nucleon potential, can be parametrized in terms of rescaling factors for the nuclear matrix 
elements, up to small residuals. Within the QRPA, the effective rescaling induced by gA variations is found to 
be almost linear, rather than quadratic in gA as naively expected. Despite this favorable feature, we find that, in 
each mechanism, the relevant lepton number violation parameter is largely degenerate with the rescaling factors; in 
particular, for gA € [0.8, 1.27] the total mpp {Mpp) uncertainty in numerical experiments amounts to a factor of about 
two (three). Moreover, the light and heavy neutrino exchange mechanisms turn out to be largely indistinguishable 
from a phenomenological viewpoint. The stringent conditions needed to lift the degeneracies between particle and 
nuclear physics aspects of Oz'/3/3 decay have been briefly discussed. Progress may be envisaged, on the one hand, 
by studying further decay mechanisms and candidate nuclei and, on the other hand, by understanding the various 
theoretical uncertainties associated to the QRPA and other nuclear models, and especially the effective functional 
dependence of the matrix elements on gA- 
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